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Abstract 

We consider the transverse-momentum (g-p) distribution of Standard Model Higgs 
bosons produced by gluon fusion in hadron collisions. At small qt {qr ^ 
nT'H, being the mass of the Higgs boson), we resum the logarithmically- 

enhanced contributions due to multiple soft-gluon emission to all order in QCD 
perturbation theory. At intermediate and large values of {It ~ "^//), we consistently 
combine resummation with the known fixed-order results. We use the most 
advanced perturbative information that is available at present: next-to-next-to- leading 
logarithmic resummation combined with the next-to-leading fixed-order calculation. 
We extend previous results including exactly all the perturbative terms up to order ag 
in our computation and, after integration over g^-, we recover the known next-to-next- 
to-leading order result for the total cross section. We present numerical results at the 
Tevatron and the LHC, together with an estimate of the corresponding uncertainties. 
Our calculation is implemented in an updated version of the numerical code HqT. 
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1 Introduction 



One of the major tasks of the physics program at high-energy hadron colhders, such as the Fermilab 
Tevatron and the CERN LHC, is the search for the Higgs boson and the study of its properties. 

Gluon-gluon fusion, through a heavy-quark (mainly top-quark) loop, is the main production 
mechanism of the Standard Model (SM) Higgs boson over the entire range of Higgs boson masses 
(100 GeV ^ruH ^ 1 TeV) to be investigated at the LHC. At the Tevatron the gluon fusion process, 
followed by the decay H — )■ WW — )■ l^l~vv>, gives the dominant contribution to the Higgs signal 
in the range of mass 140 GeV ^itlh ^ 180 GeV. In this mass region, first constraints beyond the 
LEP lower bound of 114.4 GeV [I] were established: the SM Higgs boson was excluded at 95% 
confidence level by CDF and DO collaborations in the mass range 156 GeV < mu < 177 GeV [2]. 
The first results of the ATLAS and CMS collaborations presented at EPS 2011 conference [3], and 
updated for Lepton Photon 2011 [3], dramatically extend the excluded region over most of the 
mass range between 145 and 466 GeV. 

The above exclusion relies on accurate theoretical predictions [5, 6J for the inclusive gg ^ H 
cross section, which is now known up to next-to- next-to-leading order (NNLO) [7J, with the 
inclusion of soft-gluon contributions up to next-to-next-to-leading logarithmic accuracy (NNLL) 
[8], and two- loop electroweak effects [9] 0. 

In this paper we consider the transverse momentum [qt) spectrum of the SM Higgs boson 
H produced by the gluon fusion mechanism. This observable is of direct importance in the 
experimental search. A good knowledge of the qx spectrum can help to set up strategies to 
improve the statistical significance. When studying the qt distribution of the Higgs boson in 
QCD perturbation theory it is convenient to define two different regions of qt- In the large-g^ 
region (g^ ~ ttlh), where the transverse momentum is of the order of the Higgs boson mass ttlh, 
perturbative QCD calculations based on the truncation of the perturbative series at a fixed order 
in as are theoretically justified. In this region, the QCD radiative corrections are known up to 
the next-to-leading order (NLO) [HI |T2l [13] and QCD corrections beyond the NLO are evaluated 
in Ref. pJ:J, by implementing threshold resummation at the next-to-leading logarithmic (NLL) 
level. 

In the small-gT region {qt ^ Tan), where the bulk of the events is produced, the convergence 
of the fixed-order expansion is spoiled by the presence of large logarithmic terms, a'^\n"^{m'jj/q^). 
To obtain reliable predictions, these logarithmically-enhanced terms have to be systematically 
resummed to all perturbative orders [ISl [ISl [IB dSl [TS]. It is then important to consistently 
match the resummed and fixed-order calculations at intermediate values of qt, in order to obtain 
accurate QCD predictions for the entire range of transverse momenta. 

The resummation of the logarithmically enhanced terms is effectively (approximately) per- 
formed by standard Monte Carlo event generators. In particular, MC@NLO [20] and POWEG 
[2T] combine soft-gluon resummation through the parton shower with the leading order (LO) result 
valid at large qt, thus achieving a result with formal NLO accuracy. 

The numerical program HqT [TS] implements soft-gluon resummation up to NNLL accuracy 
[22] combined with fixed-order perturbation theory up to NLO in the large-gr region [T3|. The 

'''Updated predictions for the inclusive Higgs production cross sections at the LHC are presented in Ref. |10) . 
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program is used by the Tevatron and LHC experimental collaborations to reweiglit the spectrum 
of the Monte Carlo event generators used in the analysis and is thus of direct relevance in the 
Higgs boson search. 

The program HqT is based on the transverse-momentum resummation formalism described in 
Refs. [TTl [T8l [19] , which is valid for a generic process in which a high-mass system of non strongly- 
interacting particles is produced in hadron-hadron collisions. The method has so far been applied 
to the production of the SM Higgs boson [181 ISSl |2l], single vector bosons [251 126], WW [27] and 
ZZ [28] pairs, slepton pairs [29], and Drell-Yan lepton pairs in polarized collisions [30] . 

In this paper we update and extend the phenomenological analysis presented in Ref. [18]. In 
particular, we implement the exact value of the NNLO hard-coUinear coefficients l-Lf^ computed 
in Ref. [311 132], and the recently derived value of the NNLL coefficient A'^'^'^ [33] . 

We use the most advanced perturbative information that is available at present: NNLL re- 
summation at small qt and the fixed-order NLO calculation at large qt- We present numerical 
results for Higgs production at the Tevatron Run II and at the LHC and we perform a detailed 
study of the perturbative uncertainties. We also consider the normalized qx spectrum and discuss 
its theoretical uncertainties. Our calculation for the qt spectrum is implemented in the updated 
version of the numerical code HqT, which can be downloaded from [34j- Other phenomenological 
studies of the Higgs boson qr distribution, which combine resummed and fixed-order perturbative 
results at various levels of theoretical accuracy, can be found in Refs. [35] [39]. 

The paper is organized as follows. In Sect. [2] we briefiy review the resummation formalism of 
Refs. [TTl UHl nil and its application to Higgs boson production. In Sect. [3] we present numerical 
results for Higgs boson production at the Tevatron and the LHC. In Sect. IHwe summarize our 
results. 



2 Transverse- momentum resummation 



In this section we briefiy recall the main points of the transverse-momentum resummation approach 
proposed in Refs. [TTl 111 IIH]- We consider the specific case of a Higgs boson H produced by 
gluon fusion. As recently pointed out in Ref. [T9|, the gluon fusion g-r-resummation formula has 
a different structure than the resummation formula for qq annihilation. The difference originates 
from the collinear correlations that are a specific feature of the perturbative evolution of colliding 
hadron into gluon partonic initial states. These gluon collinear correlations produce, in the small- 
er region, coherent spin correlations between the helicity states of the initial-state gluons and 
definite azimuthal-angle correlations between the final-states particles of the observed high-mass 
system. Both these kinds of correlations have no analogue for qq annihilation processes in the 
small-gT region. In the case of Higgs boson production, being H a spin-0 scalar particle, the 
azimuthal correlations vanishes and only gluon spin correlations are present [T9] . 

We consider the inclusive hard-scattering process 

hi{pi) + h2{p2) H{mH,qT)+X, (1) 

where hi and h2 are the colliding hadrons with momenta pi and p2, mn and are the Higgs 
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boson mass and transverse momentum respectively, and X is an arbitrary and undetected final 
state. 



According to the QCD factorization theorem the corresponding transverse-momentum differ- 
ential cross section dau/dq^ can be written as 

aqrp Jq Jq aqj, 

(2) 

where fa/h{x, (a = q, q, g) are the parton densities of the colliding hadron h at the factorization 
scale /iF, daH,ab/dq^ are the perturbative QCD partonic cross sections, s (s = xiX2s) is the square 
of the hadronic (partonic) centre-of-mass energy, and is the renormalization scale0. 

In the region where gr ~ the QCD perturbative series is controlled by a small expansion 
parameter, as(wiH), and fixed-order calculations are theoretically justified. In this region, the 
QCD radiative corrections are known up to NLO [HI [121 US] . In the small-gr region (g^ -C mH), 
the convergence of the fixed-order perturbative expansion is spoiled by the presence of powers of 
large logarithmic terms, a^\vL^ {m\ / q^) (with 1 < m <2n — 1). To obtain reliable predictions 
these terms have to be resummed to all orders. 

We perform the resummation at the level of the partonic cross section, which is decomposed as 

,^(res.) 7^ (fin.) 
daH,ab _ (^^H,ab ^ (^^H,ab 

dq^ dq?j. dq^ 

The first term on the right-hand side contains all the logarithmically-enhanced contributions, at 
small qT, and has to be evaluated to all orders in as- The second term is free of such contributions 
and can thus be computed at fixed order in perturbation theory. To correctly take into account 
the kinematic constraints of transverse-momentum conservation, the resummation procedure has 
to be carried out in the impact parameter space b. Using the Bessel transformation between the 
conjugate variables qx and b, the resummed component da'f^J. can be expressed as 



r ^ (res.) 



dq^ 
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l-OO ^ 

{qT,mHj]asil4i),lJ.R,lJ.l) = J db - Joibqr) Wacib,mHj; as{lJ.R), fiR, fil) , (4) 



where Jo{x) is the Oth-order Bessel function. The resummation structure of W^c can be organized 
in exponential form considering the Mellin A^-moments yV§ of with respect to the variable 
z = m^j/s at fixed mH%, 

X exp{gN{as{fJ,R),L;mjj/fi'^ji,mH/Q^)} , (5) 

were we have defined the logarithmic expansion parameter L = ln[Q'^b'^ /b^), and 6o = 2e~^'^ 
ilE = 0.5772... is the Euler number). 



■I-Throughout the paper we use parton densities f{x,fi^) and running coupling as{fJ.%) as defined in the MS 
scheme. 

§For the sake of simplicity we are presenting the resummation formulae only for the specific case of the diagonal 
terms in the flavour space. In general, the exponential is replaced by an exponential matrix with respect to the 
partonic indeces (a detailed discussion of the general case can be found in Ref. 18 ). 
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The scale Q ~ mn, appearing in the right-hand side of Eq. ([S]), named resummation scale 
parameterizes the arbitrariness in the resummation procedure. As a matter of fact the argument 
of the resummed logarithms can always be rescaled as ln(m|^6^) = ln((5^6^) + ln(m|^/(5^) (as long 
as Q ~ niH and independent of b). Although (i.e., the product Ti^ x expj^Ar}) does not 
depend on Q when evaluated to all perturbative orders, its explicit dependence on Q appears when 
W§ is computed by truncation of the resummed expression at some level of logarithmic accuracy 
(see Eq. below). As in the case of fin and fip, variations of Q around rriH can thus be used to 
estimate the uncertainty from yet uncalculated logarithmic corrections at higher orders. 



The form factor expj^jv} is universal (process independent)0 and contains all the terms OgL™ 
with 1 < m < 2n, that order-by-order in as are logarithmically divergent as 6 — >■ oo (or, equiva- 
lently, g-r — ?■ 0). Furthermore, due to the exponentiation piopertj, all the logarithmic contributions 
to Qn with n + 2 < m < 2n are vanishing. The exponent Qn can be systematically expanded as 

gN{as,L;mjj/iJ%mjj/Q^) = L g^-^^asL) + g^^\asL] m]j / m]j / Q'^) 

+ '^gfia^L- ml/fi% m\IQ^) + 0{alL^~^) (6) 

TT 

where the term Lg^^^ resums the leading logarithmic (LL) contributions agL"+^, the function g^^ 
includes the NLL contributions OgL" g^^ controls the NNLL terms d^LP'~^ |22l [33] and so 
forth. The explicit form of the functions g'^'^\ g^'^ and g^^ can be found in Ref. [T8] . 

The process dependent function 'K^ does not depend on the impact parameter h and it includes 
all the perturbative terms that behave as constants as 6 — )■ oo. It can thus be expanded in powers 
of as = asif^n)- 

L vr 

+ K'^"\<l^^l,<l^^l.<IQ') + o{al^ , (7) 

where a^^P {as,mH) is the partonic cross section at the Born level. The first order Ti^'^^'' [IQ] and 
the second order 'H^'^'^^ [311 [32] coefficients in Eq. ([7]), for the case of Higgs boson production in 
the large- approximation, are known. 

To reduce the impact of unjustified higher-order contributions in the large-g^ region, the loga- 
rithmic variable L in Eq. which diverges for 6 — )■ 0, is actually replaced by L = In {Q'^b'^/bl + 1) 
[T8| [23] . The variables L and L are equivalent when Qb ^ 1 (i.e. at small values qt), but they 
lead to a different behaviour of the form factor at small values of b. An additional and relevant 
consequence of this replacement is that, after inclusion of the finite component (see Eq. ([H])), we 
exactly recover the fixed-order perturbative value of the total cross section upon integration of the 
qt distribution over qt (i.e., the contribution of the resummed terms vanishes upon integration 
over qt). 

The finite component of the transverse-momentum cross section da^^^'^ (see Eq. (j3])) does not 
contain large logarithmic terms in the small-gr region, it can thus be evaluated by truncation of 



^It only depends on the partonic channel that produces the Born cross section. It is thus usually called quark 
or gluon Sudakov form factor. 
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the perturbative series at a given fixed order. In practice it is computed as follows 



-da 



(fin.) . 
H,ab 



-da 



dq, 



f.o. 



dq, 



H,ab 
T - 



-da 



(res.) 
H,ab ' 



f.o. 



dq. 



f.o. 



(8) 



where we have introduced the subscript f.o. to denote the perturbative truncation of the various 
terms. This matching procedure combines the resummed and the finite component of the partonic 
cross section by avoiding double-counting in the intermediate gr-region and allows us to achieve 
a prediction with uniform theoretical accuracy over the entire range of transverse momenta. 

In summary, to carry out the resummation at NLL+LO accuracy, we need the inclusion of 
the functions g''^\ g^^\ 'H^'^^\ in Eqs. f l6|7|) . together with the evaluation of the finite component 
at LO (i.e. at 0{as)) in Eq. ([8]); the addition of the functions g]^ and Hj^ , together with the 
finite component at NLO (i.e. at C(ag)) leads to the NNLL+NLO accuracy^. We point out that 
our best theoretical prediction (NNLL+NLO) includes the full NNLO perturbative contribution 
in the small-gr region plus the NLO correction at large-gy- In particular, the NNLO result for the 
total cross section is exactly recovered upon integration over qx of the differential cross section 
dau/dqT at NNLL+NLO accuracy. 

Finally we recall that the resummed form factor exp{^7v(«s(/^|;)) L)} has a singular behaviour, 
related to the presence of the Landau pole in the QCD running coupling, at the values of h where 
as{iJ^%)L > 7r//3o (/3o is the first-order coefficient of the QCD 13 function). To perform the inverse 
Bessel transformation with respect to the impact parameter b a prescription is thus necessary. We 
deal with this singularity by using the regularization prescription of Refs. |l2l HS] : the singularity 
is avoided by deforming the integration contour in the complex b space. 



3 The qt spectrum of the Higgs boson at the Tevatron and 
the LHC 

In this section we consider Higgs boson production by gluon fusion at the Tevatron {^/s = 1.96 
TeV) and the LHC (a/s = 7 TeV and 14 TeV). We present our resummed results at NNLL+NLO 
accuracy, and we compare them with the NLL+LO results. For the Tevatron we choose ttlh = 165 
GeV. For the LHC at = 7 and = 14 TeV we fix mn = 165 GeV and m// = 125 GeV, 
respectively. 

The results we present in this section are obtained with an updated version of the numerical 
code HqT [M]. The new version of this code was improved with respect to the one used in Ref. |18j . 
The main differences regard the implementation of the second-order coefficients T-Lj^ computed 
in Ref. [31] (the numerical results in Ref. [18] were obtained by using a reasonable approximation 
of this coefficient) and the use of the recently derived value of the coefficient A^^^ |33j which 
contributes to the NNLL function g^^^ (the results in Ref. [18] were obtained by using the A^^'^ 
value from threshold resummation [S]). We have checked the quantitative effect of the exact 
values of 'H^'^'^^ and A^^^ at the Tevatron and the LHC. We find that the effect is generally small 

"The evaluation of the second-order coefficient H^'^'^^ for complex values of iV, necessary to perform the inverse 
Mellin transform, is obtained using the numerical results of Ref. [H] . 
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(at the level of about 1 — 2% at the LHC at 14 TeV, 2 — 3% at the Tevatron, and at the LHC with 
7 TeV). We also find that the exact values of H^'*-^-* and A^^^ have the same qualitative impact: 
it makes the g^-spectrum (slightly) harder. 

The calculation is performed strictly in the large- approximation. 

The hadronic qt cross section at NNLL+NLO (NLL+LO) accuracy is computed by using 
NNLO (NLO) parton distributions functions (PDFs) with as(yu|.) evaluated at 3- loop (2-loop) 
order. This choice of the order of the parton densities and as is fully justified both in the small-g^ 
region (where the calculation of the partonic cross section includes the complete NNLO (NLO) 
result and is controlled by NNLL (NLL) resummation) and in the intermediate-gr region (where 
the calculation is constrained by the value of the NNLO (NLO) total cross section). Recent sets of 
parton densities, which are obtained by analyses of various collaborations, are presented in Refs. 
|1S1 SSI SZl SSI Sn]- Since the main purpose of our work is the study of the qt distribution up 
to the NNLL+NLO, we consider here only the PDFs sets of Refs. [ll, SH SSI SS] , which provide 
NNLO parton densities with Nf = 5 (effectively) massless quarks. Moreover, to avoid multiple 
presentations of similar results, we use the MSTW2008 parton densities unless otherwise stated 
(the results in Ref. [18] were obtained by using the MRST2004 set [50J). 

As discussed in Sect. [2], the resummed calculation depends on the factorization and renormal- 
ization scales and on the resummation scale Q. Our convention to compute factorization and renor- 
malization scale uncertainties is to consider independent variations oi hf and fiR by a factor of two 
around the central values /iF = yU/j = (i-e. we consider the range mn/^ < {/ii?,/^/?} < 2m//), 
with the constraint 0.5 < fJ^p/f^R < 2. Similarly, we follow Ref. ^18j and choose Q = mn/^ as 
central value of the resummation scale, considering scale variations in the range mn/^ < Q < rriH- 

In Fig. [U (left panels) we present the NLL+LO qt spectrum of a Higgs boson at the Tevatron, 
and at the LHC with ^/s = 7 TeV and ^/s = U TeV. The NLL+LO result (solid lines) at 
the default scales {fip = fJ'R = Q = rnn/^) are compared with the corresponding LO results 
(dashed lines). The LO finite component of the spectrum (see Eq. ([3])) is also shown for comparison 
(dotted lines). We see that the LO result diverges to +oo as qt — ?■ 0. The resummation of the 
small-gr logarithms leads to a well-behaved distribution: it vanishes as qx — ?■ 0, has a kinematical 
peak, and tends to the corresponding LO result at large values of qt- The finite component 
smoothly vanishes as qt ^ but gives a sizable contribution to the NLL+LO result in the low-gr 
region. 

The results in the right panels of Fig. [1] are analogous to those in the left panels although sys- 
tematically at one order higher. The qt spectrum at NNLL+NLO accuracy (solid line) is compared 
with the NLO result (dashed line) and with the NLO finite component of the spectrum (dotted 
line). The NLO result diverges to —oo as qx ^ and, at small values of qt, it has an unphysical 
peak (the top of the peak is above the vertical scale of the plot) that is produced by the numerical 
compensation of negative leading and positive subleading logarithmic contributions. In the region 
of intermediate values of qt (say, around 50 GeV), the difference between the NNLL+NLO and 
NLO results gives a sizable contribution with respect to the NLO finite component. This differ- 
ence is produced by the logarithmic terms (at NNLO and beyond NNLO) that are included in 
the resummed calculation at NNLL accuracy. At large values of qt the contribution of the NLO 
finite component noticeably increases. This behaviour indicates that the logarithmic terms are no 
longer dominant and that the resummed calculation cannot improve upon the predictivity of the 
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Figure 1: The qt spectrum of Higgs bosons at the Tevatron and the LHC. Results shown are at 
NLL+LO (left panels) and NNLL+NLO (right panels) accuracy. Each result is compared to the 
corresponding fixed- order result (dashed line) and to the finite component (dotted line) in Eq. 
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fixed-order expansion. 



Comparing the left and riglit panels of Fig. [H we see that the size of the qt spectrum increases 
at NNLL+NLO accuracy with respect to the NLL+LO accuracy. The height of the peak at 
NNLL+NLO is larger than at NLL+LO. The NNLO total cross section, which fixes the value of 
the qt integral of our NNLL+NLO result, is larger than the NLO total cross section (by about 
30% at the Tevatron and 25% at the LHC). This is due to the positive contribution of both the 
NNLO terms at small qt (the "H^'^^^ coefficient of the the "H^ function and the g^^^ function in 
the Sudakov form factor) and the NLO finite component at intermediate and large values of qt- 



Comparing Fig. l(a)[l(b) with Fig. 1(c), 1(d) and Fig 1(e), 1(f) we see that the spectrum is 



harder at the LHC than at the Tevatron. The peak of the NNLL+NLO curve moves from qx 
GeV at the Tevatron, to qt ~ 10 GeV at the LHC at y/s = 7 TeV, to qr ~ 12 GeV at the LHC 
at = 14 TeV. 

In Fig. [2]we show the scale dependence of the NLL+LO (dashed lines) and NNLL+NLO (solid 
lines) results. In the left panels we consider variations of the renormalization and factorization 
scales. The bands are obtained by varying fi^ and fip as previously described in this section. 
We note that, in the region of small and intermediate transverse momenta (gT~70 GeV), the 
NNLL+NLO and NLL+LO bands overlap. This feature, which is not present in the case of 
the fixed-order perturbative results at LO and NLO, confirms the importance of resummation to 
achieve a stable perturbative prediction. In the region of small and intermediate values of qt, 
we observe a sensible reduction of the scale dependence going from NLL+LO to NNLL+NLO 
accuracy. At the peak the reduction is from ±20% to ±13% at the Tevatron, and from ±11% to 
±8% (±12% to ±7%) at the LHC with = 7 (^i = 14) TeV. Although /^r and /xp are varied 
independently, we find that the dependence on fi^ dominates at any value of qt- 

We point out that the Qt region where resummed perturbative predictions are definitely signif- 
icant is a wide region from intermediate to relatively-small (say, close to the peak of the distribu- 
tion) values of qx- In fact, at very small values of qx (e.g. qr ^ 10 GeV) the size of non-perturbative 
effects is expected to be importanlF*!. while in the high-g^ region (e.g. qr <^ rnn GeV) the resum- 
mation of the logarithmic terms cannot improve the predictivity of the fixed-order perturbative 
expansion. The inset plots in the figure show the region from intermediate to large values of g-r- 
At large qx-, the NLL+LO and NNLL+NLO results deviate from each other, and the deviation 
increases as gr increases. As previously stated, this behaviour is not particularly worrying since, 
in the large-gr region, the resummed results loose their predictivity and should be replaced by 
customary fixed-order results. 

In the right panels of Fig. [2] we consider resummation scale variations. The bands are obtained 
by fixing /i/j = = itih and varying Q between niH/A and mu- Performing variations of 
the resummation scale, we can get further insight on the size of yet uncalculated higher-order 
logarithmic contributions at small and intermediate values of gr- We find that, in the region of 
the peak, at the Tevatron the scale dependence at NNLL+NLO (NLL+LO) is about ±4% (±10%). 
At the LHC with y/s = 1 TeV the scale dependence at NNLL+NLO (NLL+LO) is about ±3% 
(±8%) and at = 14 it is about ±3% (±13%). 

Comparing the left and right panels of Fig. [2], we see that, in the small and intermediate 



*See the discussion at the end of this Section. 
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Figure 2: The Qt spectrum of Higgs bosons at the Tevatron and the LHC. The bands are obtained 
by varying and fiR (left panels) and Q (right panels) as described in the text. 
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qt region, at NNLL+NLO accuracy, the factorization and renormalization scale dependence is 
definitely larger than the resummation scale dependence. 



The integral over qt of the resummed NNLL+NLO (NLL+LO) spectrum is in agreement (for 
any values of /i/j, and Q) with the value of the corresponding NNLO (NLO) total cross section 
to better than 1%, thus checking the numerical accuracy of the code. We also note that the large- 
ly region gives a little contribution to the total cross section; therefore, the total cross section 
constraint mainly acts as a perturbative constraint on the resummed spectrum in the region from 
intermediate to small values of gr- 
in Fig. |3] (left panels) we report our NLL+LO and NNLL+NLO total scale uncertanty bands 
(the inset plots show the large-gr region). The bands represent our best estimate of the pertur- 
bative uncertainty, and they are obtained by performing scale variations as follows. We indepen- 
dently vary jip^l^R and Q in the ranges mull < {/iF^/^i?,} < 2mH and mu/A < Q < ttih, with 
the constraints 0.5 < fJ^p/f^R < 2 and 0.5 < Q/fiR < 2. The constraint on the ratio ^if/i^r is the 
same as used in Fig. |2l it has the purpose of avoiding large logarithmic contributions (powers of 
ln(/x|,//x|j)) that arise from the evolution of the parton densities. Analogously, the constraint on 
the ratio Q/ jiR avoids large logarithmic contributions (powers of ln(Q^//i^)) in the perturbative 
expansion of the resummed form factor^ expj^Ar} (see Eq. ([6])). We recall (see e.g. Eq. (19) of 
Ref. y^) that the exponent ^at of the form factor is obtained by integration of perturbative 
functions of as(g^) over the range h'^/h'^ < g^ < Q^- To perform the integration with system- 
atic logarithmic accuracy, the running coupling Q;s(g^) is then expressed in terms of as{^R) (and 
ln(g^//i|,)). As a consequence, the renormahzation scale fiR should not be too different from the 
resummation scale Q, which controls the upper bound of the integration. 

A more effective way to show the perturbative uncertainties is to consider the fractional differ- 
ence with respect to a 'reference' central prediction. We choose the NNLL+NLO result at central 
value of the scales as 'reference' result, Xq, and we show the ratio {X — Xc)/Xc in Fig. [3] (right 
panels). The label X refers to the NNLL+NLO results including scale variations (solid lines), and 
to the NLO results including scale variations (dashed lines). 

We comment on the overall perturbative uncertainty band of our results in Fig. [3] starting from 
the Tevatron. The NNLL +NLO (NLL+LO) uncertainty is about ±13% (±28%) at the peak, it 
decreases to about ±10% (±23%) in the region up to qt = 30 GeV, and becomes ±18% (±20%) 
at qr = 60 GeV. In the region beyond qr ~ 80 GeV the resummed result looses predictivity, and 
its perturbative uncertainty becomes large. 



In Fig. 3(b) the scale variation band of the NLO result is compared to the NNLL+NLO band. 
The NLO band is obtained by varying fip and fipt as for the NNLL+NLO calculation (the NLO 
calculation does not depend on the resummation scale Q). We see that at large values of qx the 
NLO and NNLL+NLO bands overlap, and the NLO result has smaller uncertainty. As qx becomes 
smaller than about 80 GeV, the NNLL+NLO has a smaller uncertainty, and the bands marginally 
overlap. In this region of transverse momenta, the effect of resummation starts to set in. When 
qT becomes smaller and smaller, the NLO band quickly deviates from the NNLL+NLO band and 
the NLO result becomes unreliable. 

We now consider the perturbative uncertainty at the LHC, a/s = 7 TeV. The NNLL +NLO 



^^We do not apply additional constraints on the ratio Q/iip, since the form factor does not depend on fj,p. 
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Figure 3: The qt spectrum of Higgs bosons at the Tevatron and the LHC: NNLL+NLO (solid) 
and NNL+LO (dashes) uncertainty bands (left panels); NNLL+NLO (solid) and NLO (dashes) 
uncertainty bands relative to the central NNLL+NLO result (right panels). 
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(NLL+LO) uncertainty is about ±10% (±22%) at the peak, it decreases to about ±8% (±19%) 
in the region up to qx = 30 GeV, and becomes ±10% (±18%) at qt = 60 GeV. In the region 
beyond qx ~ 120 GeV the resummed result looses predictivity, and its perturbative uncertainty 



becomes large. In Fig. 3(d) we compare the NLO and NNLL+NLO bands. The qualitative 



features are similar to Fig. 3(b) at large values of qt the NLO and NNLL+NLO scale uncertainty 
bands overlap, and the NLO result has smaller uncertainty. As qt becomes smaller than about 
120 GeV, the NNLL+NLO has a smaller uncertainty, but the bands still overlap. In the region 
of intermediate transverse momenta {qx ~ 50 GeV), the bands marginally overlap and the NLO 
result underestimates the cross section. When qt becomes smaller, the NLO band quickly deviates 
from the NNLL+NLO band and the NLO result becomes unreliable. 

We finally consider the perturbative uncertainty at the LHC when ^/s = 14 TeV. The NNLL 
+NLO (NLL+LO) uncertainty is about ±9% (±25%) at the peak, it decreases to about ±8% 
(±19%) in the region up to qx = 30 GeV, and moves to ±12% (±19%) at qr = 60 GeV. In 
the region beyond qr ~ 150 GeV the resummed result looses predictivity, and its perturbative 



uncertainty becomes large. In Fig. 3(f) we compare the NLO and NNLL+NLO scale uncertainty 



bands. The qualitative features are similar to those of Figs. 3(b), 3(d): at large values of qx 
the NLO and NNLL+NLO bands overlap and the NLO result has smaller uncertainty. In the 
region of intermediate transverse momenta {qx ~ 50 GeV), the bands marginally overlap and the 
NLO result underestimates the cross section. When qx becomes smaller, the NLO result becomes 
unreliable. 



Comparing Fig. [3(a)p(b)] with Fig. [3(c)p(d)] and Fig. [3(ip(f) we see that perturbative un- 



certainties are larger at the Tevatron than at the LHC. We also note that our NNLL+NLO result 
is much more stable at the LHC than at the Tevatron, where its validity is confined to a smaller 
region of transverse momenta. This is not completely unexpected. At smaller values of the center 
of mass energy, the production of the Higgs boson is accompanied by softer radiation, and thus 
the qx spectrum is softer than at the LHC. 

We conclude this section with a discussion on the uncertainties on the normalized qx spectrum 
(i.e., l/o" X da/dqx). As mentioned in the introduction, the typical procedure of the experimental 
collaborations is to use the information on the total cross section |10j to rescale the best theoretical 
predictions of Monte Carlo event generators, whereas the NNLL+NLO result of our calculation, 
obtained with the public program HqT, is used to reweight the transverse-momentum spectrum of 
the Higgs boson obtained in the simulation. Such a procedure implies that the important infor- 
mation provided by the resummed NNLL+NLO spectrum is not its integral, i.e. the total cross 
section, but its shape. The sources of uncertainties on the shape of the spectrum are essentially the 
same as for the inclusive cross section: the uncertainty from missing higher-order contributions, 
estimated through scale variations, and PDF uncertainties. One additional uncertainty in the qx 
spectrum that needs be considered comes from Non-Perturbative (NP) effects. 

We remind the reader that the quantitative predictions presented in this paper are obtained 
in a purely perturbative framework. It is known that the transverse-momentum distribution 
is affected by NP effects, which become important as qx becomes small. A customary way of 
modelling these effects is to introduce an NP transverse-momentum smearing of the distribution. 
In the case of resummed calculations in impact parameter space, the NP smearing is implemented 
by multiplying the 6-space perturbative form factor by an NP form factor. The parameters 
controlling this NP form factor are typically obtained through a comparison to data. Since there 
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is no evidence for the Higgs boson yet, the procedure to fix the NP form factor is somewhat 
arbitrary. Here we follow the procedure adopted in Ref . |T8] , and we multiply the resummed form 
factor in Eq. (H]) by a gaussian smearing S^p = exp{—gb'^}, where the parameter g is taken in the 
range {g = 1.67 — 5.64 GeV^) suggested by the study of Ref. |37j^ . The above procedure can give 
us some insight on the quantitative impact of these NP effects on the Higgs boson spectrum. 

In Fig. m (left panels) we compare the NNLL+NLO shape uncertainty as coming from scale 
variations (solid lines) to the NP effects (dashed lines). The bands are obtained by normaliz- 
ing each spectrum to unity, and computing the relative difference with respect to the central 
normalized prediction obtained with the MSTW2008 NNLO set (with g = 0). A comparison 
of Fig. 4(a)|4(c)|4(e) to Fig. 3(b)|3(d)|3(f) shows that the scale uncertainty on the normalized 
NNLL+NLO distribution is smaller than the corresponding uncertainty on the NNLL+NLO re- 
sult. This is not unexpected: a sizeable contribution to the uncertainties shown in Fig. Incomes 
actually from uncertainties on the total cross section, which do not contribute in Fig. |H In 
other words, studying uncertainties on the normalized distribution allows us to assess the true 
uncertainty in the shape of the resummed qt spectrum. 



At the Tevatron (Fig. 4(a) ) such scale uncertainty ranges from +8% — 3% in the region of the 
peak, to +3% — 8% when qt ~ 50 GeV. At larger values of qr the uncertainty of the NNLL+NLO 
resummed distribution increases consistently with the behaviour observed in Fig. 3(b) The in- 



clusion of the NP effects makes the distribution harder, the effect ranging from 10% to 20% in the 
very small-gr region. For qt ^ 10 GeV the impact of NP effects is of the order of about 5% and 
decreases as qt increases. At the LHC, y/s = 7 TeV (Fig. 4(c) ) the scale uncertainty ranges from 
+5% - 3% in the region of the peak to +5% - 4% at qt ~ 80 GeV. At the LHC, = U TeV 
(Fig. 4(e) ) the shape uncertainty ranges from +5% — 3% in the region of the peak to +8% — 9% 
at qt ~ 100 GeV. The impact of NP effects is similar at ^/s = 7 and 14 TeV: it ranges from about 
10% to 20% in the region below the peak, is about 3 — 4% for qx ~ 20 GeV, and quickly decreases 
as qt increases. We conclude that the uncertainty from unknown NP effects is smaller than the 
scale uncertainty, and is comparable to the latter only in the very small qt region. 

The impact of PDF uncertainties at 68% CL on the shape of the qt spectrum is studied in 
Figs. [4(b)|4(d)[4(f)[ B y evaluating PDF uncertainties with MSTW2008 NNLO PDFs (red band 
in Figs. 4(b)[4(d)[4(f) ) we see that the uncertainty is at the ±1 — 2% level, both at the Tevatron 



and at the LHC. The use of different PDF sets affects not only the absolute value of the NNLO 
cross section (see e.g. Ref. [51]) but also the shape of the qt spectrum. The predictions obtained 
with NNPDF 2.1 PDFs are in good agreement with those obtained with the MSTW2008 set 
and the uncertainty bands overlap over a wide range of transverse momenta. On the contrary, 
the prediction obtained with the ABKM09 NNLO set is softer and the uncertainty band does 
not overlap with the MSTW2008 band. This behaviour is not completely unexpected: when the 
Higgs boson is produced at large transverse momenta, larger values of Bjorken x are probed, where 
the ABKM gluon is smaller than MSTW2008 one. The JR09 band shows a good compatibility 
with the MSTW2008 result, at least at the Tevatron and at the LHC for ^/s = 7 TeV, where 
the uncertainty is however rather large. At the LHC for ^/s = 14 TeV the differences with the 
MSTW2008 result are more pronounced. 



WWe note that the inclusion of this smearing factor does not change the overaU normahzation, since SNp{b 
0) = 1 
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Figure 4: Uncertainties in the normalized Qt spectrum of the Higgs boson at the Tevatron and the 
LHC. Left panels: the NNLL+NLO uncertainty bands (solid) computed as in Fig. 3 compared to 
an estimate of NP effects (dashed). Right panels: PDF uncertainties bands at 68% CL. All results 
are relative to the NNLL+NLO central value computed with MSTW2008 NNLO PDFs. 
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4 Summary 



In this paper we have considered the qt spectrum of Higgs bosons produced in hadron coUisions, 
and we have presented a perturbative QCD study based on transverse-momentum resummation 
up to NNLL+NLO accuracy. 

We have followed the formalism developed in Refs. [171 [TBI [I9] , which is valid for the production 
of a generic high-mass system of non strongly-interacting particles in hadron collisions. The 
formalism combines small-g^ resummation at a given logarithmic accuracy with the fixed-order 
calculations. It implements a unitarity constraint that guarantees that the integral over of the 
differential cross section coincides with the total cross section at the corresponding fixed-order 
accuracy. This leads to QCD predictions with a controllable and uniform perturbative accuracy 
over the region from small up to large values of qt- At large values of qt, the resummation 
formalism is superseded by customary fixed-order calculations. 

We have considered Higgs bosons produced by gluon fusion in pp collisions at the Tevatron and 
pp collisions at LHC energies, and we have presented an update of the phenomenological analysis 
of Ref. [18]. The calculation now includes the exact value of the NNLO hard-collinear coefficients 
y.^^'^^ computed in Ref. [SUES], and the recently derived value of the NNLL coefficient A^^'^ [55] . 

We have performed a study of the scale dependence of our results to estimate the corresponding 
perturbative uncertainty. In a wide region of transverse momenta the size of the scale uncertainties 
is considerably reduced in going from NLL+LO to NNLL+NLO accuracy. 

Our calculation for the qt spectrum is implemented in the updated version of the numerical 
code HqT. We have argued that, given the use that is currently done of our numerical program, the 
important information is in the shape of the qx spectrum. We have thus studied the uncertainties 
of the normalized spectrum, comparing scale and PDF uncertainties, and estimating the impact 
of NP effects. 
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